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ABSTRACT 

Acetic acid (CH3COOH) has been detected mainly in hot molecular cores where the 
distribution between oxygen (O) and nitrogen (N) containing molecular species is co- 
spatial within the telescope beam. Previous work has presumed that similar cores with 
co-spatial O and N species may be an indicator for detecting acetic acid. However, does 
this presumption hold as higher spatial resolution observations become available of large 
O and N-containing molecules? As the number of detected acetic acid sources is still 
low, more observations are needed to support this postulate. In this paper, we report the 
first acetic acid survey conducted with the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA) at 3 mm wavelengths towards G19.61-0.23, G29.96-0.02 and 
IRAS 16293-2422. We have successfully detected CH3COOH via two transitions toward 
G19.61-0.23 and tentatively confirmed the detection toward IRAS 16293-2422 A. The 
determined column density of CH3COOH is 2.0(1.0) x 10^^ cm^^ and the abundance ra- 
tio of CH3COOH to methyl formate (HCOOCH3) is 2.2(0.1) x 10"^ toward G19.61-0.23. 
Toward IRAS 16293 A, the determined column density of CH3COOH is ~1.6x IG^^ cm-^ 
and the abundance ratio of CH3COOH to methyl formate (HCOOCH3) is ~1.0xlO~^ 
both of which are consistent with abundance ratios determined toward other hot cores. 
Finally, we model all known line emission in our passband to determine physical condi- 
tions in the regions and introduce a new metric to better reveal weak spectral features 
that are blended with stronger lines or that may be near the l-2a detection limit. 

Subject headings: ISM: abundances — ISM: individual (G19.61-0.23, G29.96-0.02, IRAS 
16293-2422) — ISM: molecules — ISM: clouds 
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INTRODUCTION 



With high resolution obse rvations, Sgr B2 ha s been resolved into many regions of compact 
molecular /continuum emissi on (iGaume et al.lll995l): one of these regions, Sgr B2(N-LMH) (i.e., the 



"Large Molecule Heimat") (jSnvder et al 



1994 



Miao et al.lll995l ). is a molecule-rich hot core that 



has historically been the best place to search for large and complex molecules. Interstellar acetic 
acid (CH3COOH) was first detected in the hot molecular core (HMC) of Sgr B2(N-LMH) with 
the Berkeley-Illinois-Maryland Association (BIMA) and Owens Valley Radio Observatory (OVRO) 
millimeter- wavelength arrays bv lMehringer et al.l (119971). The compactness of the CH3COOH emis- 
sion in Sgr B2(N-LMH) (<3") (e.g.. iRemijan et ah 2002 ) highlights the importance of interferom- 
eters. To clearly detect the compact emission features from CH3COOH and other large molecules 
toward hot cores embedded in regions of more extended molecular emission, high resolution in- 
terferometric observations are necessar y. Thus, follow i ng th e first CH3COOH detections in Sgr 
B2(N-LMH) and using the BIMA array, i Remiian et aD (|2002l ) confirmed CH3COOH in Sgr B2(N- 
LMH) and discovered the second CH3COOH source, the hot core toward W51e2. 



Remijan et al.l (|2003l ) launched an extensive CH3COOH survey of 12 galactic HMCs, including 
high-mass (> 10 Mq) and low-mass (< 10 Mq) sources, and they found the thir d CH3COOH 



source toward the high-mass HMC G34.3+0.15. Using the IRAM 30-m telescope, ICazaux et al 
(I2OO3I ) reported the first detection of CH3COOH in IRAS 16293-2422 via the 9*,9 - 8*,8 E hue 



at 100.855 GHz. This is the first low-mass star forming region with detected CH3COOH to date, 
making only a total of four known interstellar sources of CH3COOH. At the conclusion of the 
search for new sources of CH3COOH, it was found that the three high mass regions containing 
CH3COOH emission share certain chemical properties. For example, the abundances of CH3COOH 
and its isomer methyl formate (HCOOCH3) are very differe nt; in Sgr B2(N -LMH), the abundance 
of HCOOCH3 is 26 times higher than that of CH3COOH ()Snydeii l2006l ) assuming the emission 
from CH3CO OH and HCOOCH3 a re co-spatial within the telescope beam. In addition, it has been 
suggested by lRemijan et al.l (j2004l ) that there is a relationship between the CH3COOH sources and 
the spatial separation of complex oxygen (O) (e.g. HCOOCH3 and CH3COCH3) and nitrogen (N) 
containing species (e.g. CH3CH2CN and NH2CH2CN). At the limit of the current observational 
resolution of these surveys, CH3COOH emission is only detected where the emission from O and 
N-containing species is co-spatial. In contrast, there have not been any CH3COOH detections in 
HMCs with report ed spatial separatio n of complex O a,nd N- containing species, e.g. Orion KL (e.g. 



Wyrowski et al.lll999l ). even though the observed separation 



Sutton et al.lll99,5l ) and W3(0H) (e.g. 
is very dependent on source distance. 

Gas phase and grain surface chemi stry may both play critical roles in the formation of CH3COOH 



and o ther large molecules in the ISM (jHasegawa et al 



1992 



Garrod &: Herbst 



20061 : 



Garrod et al 



20081 ). However, the exact pathways of CH3COOH formation in gas phase or on grain surfaces in 
astronomical environments are still not clear. The apparent spatial coincidence of the emission of 
complex O and N-containing spe cies with the CH3CO OH emission may suggest a strong correlation 
with CH3COOH formation (e.g.. iRemijan et al.ll2004l ). Another possibility is that CH3COOH and 
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large N-bearing species share approximately the same time scale to form in hot cores. Nevertheless, 
it is important to investigate the spatial relationship between CH3COOH and the distribution of 
other complex molecular species. Thus, we also observed methyl formate (HCOOCH3) and ethyl 
cyanide (CH3CH2CN) as proxies to determine the spatial distribution of O and N-bearing species 
toward our sources in this survey. 

To search for regions of CH3COOH emission and continue to determine its spatial relationship 
with the distribution of other complex molecular species, we have used the high sensitivity of the 
Combined Array for Research in Millimeter- wave Astronomy (CARMA) to carry out a CH3COOH 
survey toward three more hot cores. The observed sources are G19.61-0.23, G29.96-0.02 and IRAS 
16293-2422 and are described in more detail in §2. A complete description of the CARMA ob- 
servations are described in §3; the results are presented in §4; and in §5 we describe the analysis 
routine and discuss the implication of these observations with future surveys and chemical formation 
models. 



OBSERVED SOURCES 



2.1. G19.61-0.23 & G29.96-0.02 



G19.61-0.23 (hereafter, G19) is a molecule-rich ultracompact H II region, whic h includes a high- 



mass HMC. Numerous molecu les have been detected i n this regi on, including OH ([Matthews et al 



1977 : lGarav et al. Ill985ll. 



19991^. NH-^ (ICarav et al 



H9Q jCenze 



HCN (|Wu fc Evand 120031 ) and CO (|Hofner et al 



19981 ). CS (IShirlev et al 



k Downe£lll977l'l. CH-.OH dKalenskii et al.lll994l: 



Larionov et al 



20031 : IWu k Evansll2003l : iLarionov et al.lll999l 1. 
200d ). G29.96-0.02 (hereafter, G29) is another 



high-mass ultracompact region that also has been observed in many molecular species including 
CH3CN, C^^O, CH3OH, C H3CH9CN, C^-^S, CH3OCH3, HC OOC H3. ^^SO. SO^, HC 3N, H2CS, 
C2H5OH and SiO (see e.g. ICesaroni et aDll998l : loimi et al.l[2003l : beuther et ahlboO?!. and refer- 



ences therein). One of the more recent surveys toward both sources was conducted by lFontani et al 
(j2007l ) with the Institute de Radioastronomie Millimetrique 30-m (IRAM) telescope. These ob- 
servations targeted several emission features between 1 and 3 mm wavelengths of CH3CH2CN, 
vinyl cyanide (CHo CHCN) and dimethyl ether (CII3OCII3). Specifically, the observations by 
Fontani et al.l (|2007l ) were conducted to further investigate whether the so called "chemical dif- 
ferentiation" seen toward Orion and W3 is a general characteristic of high mass HMCs. However, 
no "chemical differentiation" was observed within the spatial limits defined by the size of the 
synthesized beam of their observations. 



Toward G19, from the observations of a series of CH3CH2CN transitions. iFontani et al.l (|2007l ) 
determined a rotational temperature of T,.ot=116(12) K, and from that a source averaged col- 
umn density of A^5'=2.2(0.3) x 10^^ cm^^. From a series of CH3OCH3 transitions that would be 
used to repre s ent th e excitation and distribution of other 0-bearing molecules like HCOOCH3, 
Fontani et all ((2003) determined a Trot=158(17) K, and ^¥^=2.0(0. 2) x 10^^ cm'^. Note that the 
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half power beam width (HPBW) of the IRAM 30-m telescope ranged in these 1 mm to 3 mm ob- 
servations from ~12"to ^ 22" , r espectively. In contrast, using archival data from the Submillimeter 



Array (SMA), IWu et al.l (|2009l ) found a much higher temperature for the hot core regions. Using 
methyl cyanide (CH3CN) at a spatial resolution of 2'.'6x2'.'l, they found Trot~552 K and the beam 
averaged column density, A^5'~3.4xl0^^ cm~^ without optical depth corrections. 



In a dedicated CH3COOH survey of high mass HMCs bv iRemiian et al.1 (1200J), HCOOCH3 
(methyl formate), CII3CII2CN (ethyl cyanide), CII3OII (methanol) and HCOOH (formic acid) 
were all detected at a spatial resolution of ~5" toward G19. There was also an indication of weak 
CII3COOH emission at 111.507 GHz, at the la detection limit of the BIMA array. Based on the 
detected CH3COOII sources and the upper limits from the non-detections, the abundance ratio of 
CH3COOH to HCOOCH3 seems universal in high-mass HMCs. If G19 indeed contains CH3COOH, 
the HCOOCH3 detection indicated that the CH3COOH transitions were just slightly lower than the 
BIMA sensitivity. To clearly detect the CH3COOH transition lines, we needed a higher sensitivity 
by a factor of 4, which is provided by CARMA in this survey. 



Toward G29, the observations by iFontani et al.l (|2007l ) of CH3CH2CN transitions yielded a 
Trot=121(17) K, and from that a A^5=1.5(0.3)xl0i6 cm'^. Similarly, observations of CH3OCH3 
transitions yielded a Trnj:=141(26 ) K, a nd A^s'=6.2(1.0) x 10^^ cm~^. Another recent campaign 



toward G29 was by iBeuther et al.l (j2007l ) who mapped the spatial distribution of several complex 
molecules with the SMA at a spatial resolution of ~0'.'4x0'.'3. However, no temperature or column 
density determinations were made from the emission features of these complex molecules. The only 
reported temperature of the compact emission regions was made from observations of CH3OH in 
its vibrational l y exci ted ft=l state. From the relative intensities of the detected emission features, 
Beuther et al.l (|2007l ) reported an excitation temperature of Tex~340 K with a lower limit of 220 K. 
In addition, this mapping campaign showed the distributions of CH3OH, CH3CH2CN, CH3OCH3 
and HCOOCH3 were co-spatial suggesting that G29 may be another source of CH3COOH and that 
an excitation temperature as high as 220 K may be used in the determination of the column density 
of these complex molecules (see §3). 



2.2. IRAS 16293-2422 



I RAS 16293-2422 is one of the most well-studied and observed low-mass star forming region 



(e.g., iLooney et al.l l200Cl ). It consists of two clumps denoted as component A a nd B (hereafter 



sources Aa and Ab fiChandler et al 



with two components (jLoinard et al 



I1629 3A and B). The structure of I16293A is much more com plex thari I1629 3B (e.g. I Chandler et al 
2OO5I ). I16293A cont ains two centirneter s ources Al and A2 (|Woottenlll989l ). and two submillimeter 



2OO5I). Moreover, I16293A2 may actually be a bipolar ejection 
2007. ) . Despite their low mas ses, both componen t A and B have 



beeri suggested to be as molecule-rich as high-mass HMCs (e.g. iSchoier et al.l 120021 : iCazaux et al 
2OO3I ). 
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Large 0-bearing species, including HCOOCH3 and HCOOH, and N-bearing species, including 
CH3CN (me thyl cyanide), CH2CHCN (vinyl cyanide) and CHsCHoCN, have been detected towar d 



this region (ICazaux et al.l l2003l : iBottinelli et al 



2004 : iKuan et alJ 12004 : iRemiian k HollisI I2OO6I ) . 



With the IRAM 30-m telescope, ICazaux et al.1 (|200a ) reported the first detection of CH3COOH 
in this low mass region via the 9*,9 - 8*,8 E hue at 100. 855 GHz. Thi s mak es IRAS 16293-2422 
the only low-mass CH3COOH source to date. However, ICazaux et al.l ([2003) did i iot detect the 
counterpart, 9*^9 — 8*^8 ^ line, and the 9*^9 — S,,,^ E line is known to be blended (e.g.. ICazaux et al 



2OO3I ). Comparing the IRAM single dish of a 28" beam to an estimated source size of 5" for the 



CH3COOH emission, beam dilution would severely hinder the detection of other weak CH3COOH 
lines and may even pick up the extended emission from more distributed molecular species. To 
further investigate the CH3COOH distribution toward IRAS 16293-2422, we observed two more 
CH3COOH transitions in this survey with the high resolution of CARMA, which is better coupled 
to the assumed source size and is insensitive to extended emission. 



OBSERVATIONS 



The CH3COOH survey was carried out with the CARMA A = 3 mm receivers from Spring 
2006 to Spring 2008 in the B and C configurations. The individual angular resolutions are about 
1" and 2", respectively. Table [1] summarizes the properties of the observed sources, which includes 
the source positions, flux calibrators, gain calibrators, distances, masses and LSR velocities. The 
six windows in the CARMA correlator were configured as two wide-band and four narrow-band 
windows. Each wide-band window has 500 MHz bandwidth and 15 channels while each narrow- 
band window has 31 MHz bandwidth and 63 channels. The narrow-band high spectral resolution is 
0.49 MHz per channel, which provides sufficient sensitivity to resolve the observed transitions and 
lower the line confusion effect to detect the compact emission from large molecules. The amplitude 
and phase calibrations were accomplished by observing the gain and flux calibrators. The overall 
flux calibration uncertainty is 10-15%. The CARMA flux uncertainty discussed in the text is the 
statistical uncertainty and does not account for any amplitude calibration error. In this paper, we 
will use the style of showing la statistical uncertainties in parenthesizes after the derived values. 
During data reduction, pointing and system temperatures were checked to monitor data quality. 

We observed the 10. m - 9*,9 E and A CH3COOH lines at 111.507270(20) and 111. 548533(20) 
GHz (lllvushin et al I2OO8I ) in two narrow band windows. Based on previous observations (jRemijan et al, 



2002 



200^ . they are unblended and have similar line strengths. In addition to CH3COOH, we 



included two HCOOCH3 and three CH3CH2CN lines in the other narrow-band windows. Table [2] 
summarizes the observed molecules and their transitions. All spectral line data were taken from the 
Spectral Line Atlas of Interstellar Molecules (S LAIM) (F. J. Lova s, private communication)0, the 
Cologne Database for Molecular Spectroscopy (iMiiller et al.ll2005l ) and the appropriate laboratory 



^available at |http: / /www.splatalogue.net (|Remiian et al.ll2007l l 
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references listed therein. 



4. RESULTS 



The observational data were reduced with the MIRIAD package (jSault et al.lll995l ). The final 
data presented in the figures were Hanning-smoothed over 3 channels and continuum-subtracted. 
The observations in B and C configurations were combined using natural weighting to obtain the 
best sensitivity and synthesized beams of ~2". We present the observations of each source as both 
maps of the molecular emission and the spectra at the peak of the emission (Figures 1 through 
7). The continuum emission maps shown were made from channels contained in the wideband 
windows which were free from line emission. All of the spectra from each source are shown in rows 
and columns. In each spectral figure, column (a) displays window 2 (rest frequency=107. 54375 
GHz); column (b), window 3 (rest frequency=107. 59404 GHz); column (c), window 6 (rest fre- 
quency=lll. 50727 GHz); and column (d), window 5 (rest frequency= 11 1.54853 GHz). In row 1, 
the observed spectra (Hanning smoothed over 3 channels) are overlaid on the modeled spectral line 
data represented by Gaussians (red trace). In row 2, the modeled spectral line data from row 1 
are convolved with the spectral resolution of the original observational data and co-added with 24 
mJy/beam random noise. Finally, row 3 presents the residuals from the modeled data subtracted 
from the observational data (i.e row 1 minus row 2). The maps and spectra are discussed below. 



4.1. G19. 61-0.23 

One of the most relevant discoveries in this paper is the new detection of CH3COOH toward 
G19 that highlights the importance of interferometric observations in the detection of weak spectral 
signatures confined in compact emission regions. To illustrate the distribution of the CH3COOH 
emission toward G19, we have summed the intensity of the CH3COOH lines to produce a combined 
CH3COOH contour map which is overlaid on the continuum map in grayscale in Figure [Dj^a) . Both 
the 10*, 10 — 9*^9 E and 10=,,, 10 — 9*^9 A lines were detected in the two CH3COOH spectral windows 
(Figures [2]||c) and (d)) at the canonical v^sR of 40 km s~^ for G19. While some weak unidentified 
lines may appear in the windows after Hanning smoothing, the relative line strengths of the two 
CH3COOH transitions are roughly equal as predicted by rotational spectroscopy. G19 consists 
of several high den sity components of continuum emission between 2 cm and 3 mm wavelengths 



(|Furuva et al.ll2005l ). In this case, the CH3COOH emission is located near the center of component 



C and the CH3COOH emission region has an effective scale size of 2". 

In addition to CH3COOH, we have observed transitions of HCOOCH3 and CH3CH2CN. Two 
transitions of HCOOCH3, 92,8 - 82,7 A/E and three transitions of CH3CH2CN, 124,9 - ll4,8, 124,8 - 
114,7, and 123,10 — 113,9 are shown in Figures [2|^a) and (b). All transitions are at a wl5_r=40 km 
s~^. The statistical significance of three possible unidentified lines will be discussed in Section 
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5.5. The HCOOCH3 map (Figure and CH3CH2CN map (Figure [D^c)) show emission that 

peaks primarily near component C, Hke the CH3COOH emission, yet with a south ward projection 



towa rd component A that was not seen in the lower resolution BIMA observations (jRemiian et al 



20041 ) . Table [3] summarizes the fitting results of the detected transitions towards G19. In our fitting 



routine, we assume that the rest frequencies of each molecular transition is fixed to the vlsr of the 
source and the HCOOCH3 transition that is blended with the CH3CH2CN transition shares the 
same linewidth and peak intensity as the unblended transition. Also, shown are the line widths 
and intensities of the model data for a direct comparison. 



4.2. G29. 96-0.02 

We have detected HCOOCH3 and CH3CH2CN toward G29. Figures [^a) and (b) show the 
distribution of HCOOCH3 and CH3CH2CN emission toward the G29 region, respectively. In Figure 
m we clearly detected the two blended pairs of transitions of IICOOCII3 92,8 — 82,7 E/A and 
CII3CH2CN 1249 — 114,8 and 1248 ~ 114,7- The spectrum is very similar to that of G19 (Figure 
El^a)). The CII3CH2CN 123,10 — 113,9 transition and possible U lines in Figure ID^b), also have 
similar profiles to those in Figure [2^b). 

G29 consists of two main contin uum peaks at 1.2 cni , which trace the UGH II region, and at 



862 /xm, which trace the cloud core (jBeuther et al.l 120071 ). The cloud core is the location of the 



molecular emission; however, the emission is not bright enough to contribute significantly at 3 mm 
wavelengths. In Figure [3] the continuum emission is dominated by the UG H II region. 

Figure m^c) and (d) shows the two GH3GOOH spectral windows toward G29. Unlike the data 
from G19, there are no statistically significant features at the CH3GOOH rest frequencies. There is 
a la feature that may be due to the 10=,,, 10 — 9*,9 E transition of GH3GOOH, but the corresponding 
10*, 10 — 9=K,9 A transition is noticeably absent. In the discussion, we will investigate these passbands 
more closely with respect to the modeled data. Table H] summarizes the fitting results of the 
detected transitions towards G29. The fitting criteria and model predictions are presented in the 
same manner as in Tabled 



4.3. IRAS 16293-2422 



Unlike the spectra from iGazaux et al.l (j2003l ). we resolved IRAS 16293-2422 into I16293A 
and B with our 3'.'lx l'.'9 beam (see Figure [5]). Therefore, beam dilution in this observation 
is minimized. Figure [6] and [7] shows the two GII3GOOII spectral windows toward I16293A and 
I16293B, respectively. Toward I16293A, we detect weak emission features that may be due to 
transitions of GII3GOOII but they are right at the 2a detection limit. Toward I16293B, as in 
G29, there are no statistically significant features at the GII3GOOH rest frequencies. There are la 
features that may be due to the 10,,, 10 — 9*, 9 A/E transitions in these Banning smoothed data, but 
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nothing convincing beyond that level. In §5.5 we evaluate whether the features toward I16293A 
may be due to CH3COOH. 

We have also observed two new transitions of HCOOCH3. The line profiles of HCOOCH3 shown 



in Fi gure El toward I16293A are wider and may be suffering from self absorption (jRemijan &: HoUis 



20061 ) more than the line profiles shown in Figure [7] toward I16293B, which are much more narrow. 
This suggests that I16293A and I16293B may have distinguishable molecular dynamics. Figure 
[5] (a) and (b) show the distribution of HCOOCH3 emission with respect to the 3 mm continuum 
emission toward the I16293A & B regions, respectively. Table [5] summarizes the fitting results of 
the detected transitions towards IRAS 16293-2422. The fitting criteria and model predictions are 
presented in the same manner as in Table [3l 



5. ANALYSIS & DISCUSSION 

5.1. Molecular Column Densities 

Molecular column density is one essential physical parameter needed to constrain chemical 
models (see e.g. Garrod et al. 2008 and references therein). In order to determine the column 
densities of each molecular species of interest, we assume that each region has uniform physical 
conditions, that the populations of the energy levels can be characterized by a Boltzmann distri- 
bution, and finally, that the emission is optically thin. Assuming that the molecular species is in 
local thermodynamic equilibrium (LTE) and low optical depth, the total beam averaged column 
density is 

/V rrm"^^ - 2 01 M ^n^o ^'^^^^ beam-i)^t;(km s"^) Q.e^" W/^- 
iVHcm )- 2.04x10 ^^(,rcsec2),5(GHz3) < 5^^^ Xdebye^) ' 

where is the solid angle of the beam , J lydv is the integral of the line intensity over velocity, v 
is the spectral line frequency, < Sfi"^ > is the line strength parameter, Qr is the rotational partition 
function, Tr is th e rotational temperature, and is the upper level energy of the transition 



(jMiao et al.lll995l ). In each case, we assume the source emission fills the synthesized beam of the 



observations. All of the line parameters used for the analysis are given in Table 2. 

A least-square fitting routine with Gaussian functions was used to determine the linewidths 
(Af), peak intensities (AI) and j Idv of each observed spectral line toward each source. All 
least-square fitting was done using the standard packages contained in the Mathematica software 
package. Using the spectroscopic parameters from Table 2 and by varying the temperature and 
column density over a wide range of values (e.g. the temperature range investigated was from 
10-500 K and the column density range was from 10^^ - 10^* cm~^ for each molecular species). 
Equation (1) gave a predicted integrated line intensity for each spectral line observed in our data. 
These predictions of integrated line intensity were then compared to the measured values and, 
using a least-squares fitting routine, a best fit temperature and column density was found for each 
molecule toward each source. However, there were some sources where an independent fit of the 
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temperature did not converge. These sources are described in more detail in the following sections. 
Columns 3 and 4 of Table 3, 4 and 5 give the fitted line parameters of the observed peak intensity 
and line width. Columns 5 and 6 give the modeled peak intensity based on the least square fitting 
routine and the fitted line width from the observed data, respectively. Once all the relevant line 
parameters have been determined, a model spectrum with the appropriate amount of random noise 
added and convolved with the spectral resolution of the observation can then be produced and be 
directly compared with the observed dataset (see §5.5). 

The rotational temperature generally is representative of the physical kinetic temperature 
in regions where the spatial density is larger than the critical density of the molecule in question. 
The assumption in this work is that the densities of the HMCs investigated are expected to be 
high enough for this to be true for these molecules. The rotational temperature is often found by 
observing several different transitions of a molecular species over a range of Eu values. However, 
given the limited bandwidth and low number of transitions observed in our data, it was often 
difficult to accurately const rain the temperatu re and column density of our three observed molecules 



toward these regions (e.g., ISnyder et al.ll2005l ). Therefore, when the least squares fit failed to give 
a reasonable temperature for a given molecular species, we adopted a rotational temperature from 
previous observations by other investigators toward these sources at similar spatial resolution and 
were then able to fit for the total beam averaged column densities or associated upper limits if no 
transitions were detected beyond the statistically significant la limit. It was also the case that if the 
fitting procedure was given too many free parameters, a constraint was needed on one or more or 
else it would not converge (see the G29 discussion below). To properly constrain the temperature of 



s of additiona 


transitions are necessary 


Fontani et al. 


2007 




Priedel et al. 


2004, 



and references therein). 

For the G19 region, the least square fitting routine determined a rotational temperature of 
CH3CH2CN of T^ot=161(58) K, and a beam averaged column density of iVT=6(3) x 10^^ cm'^. 
The errors on all fits are la. The routine could not find a reasonable fit for the temperature and 
column density for either HCOOCH3 or CH3COOH presumably because the detected transitions 
were too close in energy. Therefore, because the emission regions of each molecular species around 
G19 are nearly co-spatial (Figure we adopted the same rotational temperature for HCOOCH3 
and CH3COOH as found from the CH3CH2CN fit. Assuming a T^ot=161 K for HCOOCH3, we find 
a best fit to the column density of A''t=9(2) x 10^^ cm~^ and for CH3COOH, we find a best fit to the 
column density of A't'=2(1) x 10^^ cm~^. Given the minimum and maximum values of the measured 
column density, we determine an abundance ratio range of ^HCOOCHs/^CHsCOOH='i-^^- Table 
[3] gives the fitted intensities of each of the detected spectral features compared to the predicted 
intensities determined from Equation (1) for G19. 

For the G29 region, the least square fitting routine determined a rotational temperature of 
CH3CH2CN of Trot=107(49) K, and a beam averaged column density of iVT=l.l(7) x 10^^ cm~^ 



only if we adopted a rotational temperature of 150 K, determined from lOlmi et al.l (|2003l ) and 
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supported by the observations of lFontani et al.l (|2007l ) and iBeuther et al.l (|2007l ) , for the transitions 
of HCOOCH3. This temperature was also set to determine the upper hmit for the column density 
of CH3COOH. Assuming a Trot=150 K for HCOOCH3, we find a best fit to the column density 
of A^r=4(l) X 10^^ cm~^ and for CH3COOH, we find an upper limit to the column density of 
A'j'<9xl0^^ cm~^ for a relative abundance ratio upper limit of ^hcoocHs/^cHsCOOH>50- Table 
m gives the fitted intensities of each of the detected spectral features compared to the predicted 
intensities determined from equation (1) for G29. We also note the detection of another feature 
shown in Figure [5] near 107.540 GHz that is most likely due to two high energy (>110 K) spectral 
lines of CH3CH2CN. The discussion of weak spectral features near the la noise level will be given 
in section §5.5. 



A recent observation of the G19 and G29 regions with the IRAM 30-m telescope (iFontani et al 



20071 ) covered our two acetic acid line frequencies (Figure [8]). However, with a 28" beam, the lines 
detected are under the noise level of 26 mK. These lines are brighter than those observed by 
Cazaux et al.l (120031 ) . but the beam dilution is severe. The high resolution capability of CARMA 
is indeed critical for searching for acetic acid and large molecule research. 

Finally, toward IRAS 16293-2422, the least square fitting routine could not determine a ro- 
tational temperature for the detected transitions of HCOOCH3 presumably for the same reason 
as described toward G 19. Thus, we adopted a rotational temperature of 62 K determined from 
Bottinelli et al.l (120041 ). In this case, we find a best fit to the column density of HCOOCH3 of 
A'7^=1.5(3) X 10^^ cm~^ toward region A and 6(1) x 10^^ cm~^ toward region B. Since there was no 
clear detection of CH3CH2CN beyond the Icr detection limit toward either region A or B, based 
on the noise level in each of the passbands containing those lines, we determined an upper limit 
to the total beam averaged column density of CH3CH2CN to be <6xl0^^ cm^^ for region A; and 
<3xl0^^ cm~^ toward region B. For CH3COOH, the 2a detection toward region A gives a total 
beam averaged column density of ~1.6xl0^^ cm~^ and for region B, we determine and upper limit 
of <6xl0^'^ cm~^. Table [5] summarizes the fitting results of the HCOOCH3 transitions compared 
to the model intensities determined from Equation (1). 



5.2. A Proxy for the Detection of CH3COOH Sources 

We have successfully detected CH3COOH toward G19. Following Sgr B2(N-LMH), W51e2 
and G34. 3+0.15, G19.61-0.23 is the fourth high mass hot core that is a source of CH3COOH. We 
compare the CH3COOH column densities and the abundance ratios to HCOOCH3 of all detected 
CH3COOH hot cores in Table [71 The values are generally consistent. Therefore, the column 
densities and ratios can be used to constrain hot core chemical models assuming the co-spatial de- 
pendence of HCOOCH3 and CH3COOH remains. This has yet to be tested at extremely high spatial 
resolution. In addition, a consistency can be found between all the sources containing CH3COOH 
emission by comparing the relative line strengths between the detected transitions of HCOOCH3 
and CH3CH2CN to the CH3COOH emission. Taking a ratio between measured (or predicted inten- 
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sities) of the 107.594 GHz line of CH3CH2CN to the 111.507 GHz hne of CH3COOH, we find that 
this ratio is between 5.5 and 6 for those sources where CH3COOH is detected. Therefore, given 
the hne intensity of 0.21 Jy/beam for the 107.594 GHz hne of CH3CH2CN toward G29.96-0.02, the 
expected hne intensity for the 111.507 GHz hne of CH3COOH is <0.04 Jy/beam. This intensity is 
right at the la detection hmit of the current CA RMA observati o ns. T herefore, in addition to the 
criteria for the detection of CH3COOH outhned in Remijan et al. (j2003l ). we now have a good proxy 
for the detectability of the 111 GHz CH3COOH transitions relative to the detected line strengths 
of the HCOOCH3 and CH3CH2CN transitions near 107 G Hz. Further observatio ns of new sources 
of HCOOCH3 and CH3CH2CN that fit the criteria set bv iRemiian et al.l ^200^ ) and higher S/N 
observations of sources like G29, are necessary to validate this hypothesis. 



5.3. CH3COOH in IRAS 16293-2422? 

High spatial resolution maps and spectra of IRAS 16293-2422 clearly reveal that the mor- 



(Guelin et al. 


2008). 


Blake et al. 


i 



and lRodgers &: Charnleyl (120011 ) have suggested the difference in the morphologies of the emission 
of complex molecules in these regions reflects the timescales of chemical evolution. Predominantly 
N-bearing species (e.g. CH3CH2CN) represent the later stage of chemical evolution while predom- 
inantly 0-bearing species (e.g. HCOOCH3) represent the earlier stage. Their models can estimate 
the ages of hot cores by comparing this difference between the O and N morphology of emission. 
However, to date, no direct observation can accurately measure hot core ages, so the uncertainty 
of these estimates is large. 



I RAS 16293-2422 also sh ows clear distinctions in the morphology of molecular emission (IChandler et al 



2OO5I : iRemijan &: Hollisll2006l ). In fact, the relative distributions of molecular species will undoubt- 
edly give insight into their formation chemistry. It may also turn out that the relative abundances 
of molecular species and the overall population of energy levels of the detected transitions of a given 
molecule may be determined more by the formation chemistry and not on the local thermodynamic 
environment. To test this hypothesis will require high spectral and spatial resolution observations 
of a group of species that may form via a common formation pathway. 



IRAS 16293-2422 has been reported as the first low-mass CH3COOH source by lCazaux et al 



(j2003l ) using the IRAM 30-m telescope, but the one detected line was likely blended. While IRAS 
16293-2422 is not resolved in their work (beamsize of 28"), we have resolved it into I16293A and 
B with a separation of less than 5". Even though their observation used different transitions, we 
should detect our CH3COOH lines in our synthesized beam of 3'.'lxl'.'9 using their fiux estimates. 
Although at a lower fiux than expected from their results, we do detect two CH3COOH transitions 
tow ard region A but oii ly at the 2a limit. We estimate a column density of ~ 1.6 x 10^^ cm~-^. 



and Cazaux et al. (j2003l ) estimated an upper limit of 2.5 x 10^^ cm ^, which is similar. 
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Furthermore, there is no evidence beyond the la levels of any emission of CH3COOH toward 
region B. In addition, while we have observed HCOOCH3 in I16293A and B, we did not detect 
the corresponding CH3CH2CN emission features. The detection of CH3COOH toward G19 (and 
possibly toward I16293A) and the other hot cores have demonstrated that a diversity of abundant 
large O and N-bearing species are strongly correlated with CH3COOH. One possibility is that the 
mix of O and N-bearing species may serve as an efficient chemical network for CH3COOH and 
other large molecules. Thus, they need to wait for large N-bearing molecules to form first before 
they can be efficiently produced. Another possibility is that the formation of large O and N-bearing 
species is simply time dependent, and the gas-phase CH3COOH concentration just by coincidence 
starts to increase along with large N-bearing species. Given the criteria set above, it is intriguing 
that CH3COOH is detected toward I16293A given there was no clear evidence of any emission 
from CH3CH2CN. Further higher sensitivity observations are once again necessary to confirm this 
detection and determine the evolutionary stage of the I16293A and B cores. 



5.4. CH3COOH Formation 



Both gas phase and grain surface chemistry undoubtedly play important roles in hot cores (e.g. 
Hasegawa et al.lll992l : ICarrod fc: Herbstll2006l : iGarrod et al.ll2008l ). Large molecule abundances can 
be enhanced enormously by grain surface chem istry while current mode ls of gas phase chemistry 
alone are not sufficient to match observations. iGarrod &: HerbstI (|2006l ) have demonstrated that 
radicals in the icy mantles on grain surfaces can move around as the temperature gradually in- 
creases from 10 to 200 K in star forming regions. Grain surfaces become efficient environments 
for large molecules to form. However, the details of CH3COOH formation in hot cores still remain 
unclear. The correlation between CH3COOH and large N-bearing molecules may suggest that some 
N-bearing molecules act as catalysts to the CH3COOH formation. The chemical and physical evo- 
lution timescales of star formation lacks direct observational evidence to confirm chemical models. 



Without the support of observation, it is difficult to establish solid chemical models. IGarrod et al 



(|2008l ) have reported their new chemical model for CH3COOH in hot cores. Grain surface reactions 
have been enhanced to match observations. They show the time scales of hot cores warming up af- 
fect the abundances of the secondary radical CH3CO, which can react with OH to form CH3COOH. 
This effect can also cause the dissimilar abundances of structural isomers, which is consistent with 
our survey. 



5.5. Identifying Weak Spectral Features Using Model Spectra. 



When CH3COOH was first detected (jMehringer et al.l 119971 ) . it emphasized the importance 
of interferometry for the study of astrochemistry in hot cores. Star forming regions generally 
contain many regions of compact molecular emis sion. Although G1 9 has a complex structure, 
including extended and compact components (see iFuruya et al.ll2005l ). the CH3COOH detection 
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toward G19 was conducted with a high resolution beam, l'/9xl'/6, which is the smallest used 
in a CH3COOH survey to date. The three observed molecules in G19, CH3COOH, HCOOCH3 
and CH3CH2C N, are mainly locat ed at the same component, which agrees with the CH3CH2CN 
observation bv lFuruva et al.l (|2005l ). As we continue extending CH3COOH surveys to smaller hot 
cores, the required resolution of observations will be higher. 

Higher resolution observations also have the advantage of reducing the effect of line confusion 
present in a spectrum. Also, if the source size of the emitting region is well coupled to the synthe- 
sized beam of the array, weak spectral features that are lost due to beam dilution effects from single 
dish observations can be significantly resolved out of the noise. Thus, many new molecule detections 
over the last several years have used both a combina tion of single dish and array observations to 



confirm the detection (e.g.. iBelloche et al.ll2008l . l2009l ). In addition, these detections routinely use 



model spectra based on the laboratory and calculated rotational spectroscopy of energy levels and 
line strengths for molecular transitions, the physical conditions of the emitting regions including 
temperature and density and finally, the characteristics of the telescope including forward beam 
efficiency and field of view (i.e. beam size). Given the model spectra, the weak transitions of 
molecular species appear statistically significant in the observed spectrum. However, the model 
data have near infinite spectral resolution compared to the observed data and do not include an 
estimate of the noise level in the passband. These types of models are shown in the spectral line 
data presented in this work (red Gaussian traces in the first row of spectra in Figures 2, 4, 6 and 
7). 

In order to correctly ascertain whether a weak feature in the observed spectrum is due to a 
calculated line of a molecular species, we advocate taking the model data and convolving it with 
the spectral resolution of the instrument and adding the appropriate amount of random noise to the 
resultant spectrum to give a more accurate representation of the data collected by the telescope. 
The noise level is based on the measured noise of line free channels in the observed data. These 
modeled data are then subtracted from the observed data to form a residual spectrum. From these 
residuals, with the caveat of an accurate model of the physical conditions of the source, one can 
identify the remaining weak spectral features and determine if features close to the l-2a level are 
in fact statistically significant spectral line features. 

Based on the residual passbands shown in the last row of Figures 2, 4, 6 and 7, we include 
in Tables 3, 4 and 5 the possibility of unidentified lines (i.e. U-lines) in the data toward G19, 
G29 and IRAS 16293-2422. In addition, given the spectral line fit (red Gaussian traces) toward 
IRAS 16293-2422, there is a suggestion that the weak emission features shown toward I16293B 
may in fact be due to CH3COOH. However, when this model spectrum is convolved with the 
noise level in the passband, these weak features completely disappear and there are definitely no 
statistically significant feature seen in the residual passbands. As a result, while we claim that the 
emission features that are coincident with the CH3COOH transitions toward I16293A are in fact 
from CH3COOH, these features are only barely above the noise level so further high sensitivity 
observations are needed for a definitive detection. 
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It is our recommendation that if model spectra are going to be the new proxy for identifying 
weak spectral features in observational data, the model data should be first convolved to the spectral 
resolution of the telescope and the appropriate amount of noise be added to the model in order to 
make an accurate representation as possible to the observed data. Only then can the l-2cj spectral 
features be corrected assessed as statistically significant features in the observed data and not an 
artifact of the noise and spectral resolution. 



6. SUMMARY 

We have conducted a high resolution acetic acid (CH3COOH) survey at 3 mm wavelengths for 
the first time with CARMA towards two high-mass hot cores, G19.61-0.23 and G29.96-0.02, and a 
low-mass embedded protostar system, IRAS 16293-2422. We have detected CH3COOII emission in 
G19.61-0.23, via the two CII3COOII rotational transitions, 10*,io - 9*,9 E and A. While 019.61- 
0.23 consists of several clumps, high resolution mapping reveals that the OII3COOH emission is 
extremely compact (<2") toward component 0. Methyl formate (IICOOCH3) and ethyl cyanide 
(OII3OII2CN) were also detected. The emission from these two molecules is more extended than 
the OH3OOOH emission but the peaks of all three emission spectra are toward component 0. The 
CII3OOOII column density and abundance ratio with respect to methyl formate (IICOOCII3), are 
2.0(1.0) X 10^^ cm~^ and 2.2(0.1) xlO~^ respectively, which is higher but comparable to the other 
high-mass CII3OOOII sources given the current level of detection. Although acetic acid is not 
detected toward O29.96-0.02, we have detected HCOOCH3 and OH3CH2ON. The HOOOOH3 and 
CH3OH2ON emission features of O29.96-0.02 are weaker than O19.61-0.23 by a factor of 3. While 
we detect OH3COOH in O19.61-0.23, the non-detection of OH3COOH in O29.96-0.02 may be due 
to the limited sensitivity. Finally we have detected two new transitions of IICOOCII3, 92.8 ~ 82,7 E 
and A toward IRAS 16293-2422 and weak emission features at the 2a level of CII3OOOH toward 
I16293A. Th is is the first observ ation of new spectral features of OH3OOOII since the reported 



detection by lOazaux et al.l (120031 ). The HOOOOH3 spectrum toward I16293B shows narrower and 
more distinguishable line profiles than the spectrum toward I16293A, which may be suffering from 
self-absorption. However, we did not detect CII3CII2ON in IRAS 16293-2422 at the the la detection 
limit. 

The new CH3OOOH detection demonstrates the strong correlation between large O and N- 
bearing species and OH3COOII in hot cores. The timescale of hot cores may also play a role in 
0/N chemical differentiation and hence OH3COOII formation. However, we need more OH3OOOH 
detections to confirm these suggestions. The compactness of the emission from large molecules 
stresses the need of interferometers for hot core chemistry observations. Our high resolution survey 
has provided critical information about the chemical and physical properties of hot cores, which 
can be used to constrain the hot core chemical models. Higher resolution observations are still 
required to precisely determine temperature and the density distribution of large O and N-bearing 
species for physical-parameter sensitive chemical models. Through well-improved chemical models, 
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we can better understand the distinctive hot core chemistry. 
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Table 1. Observed Sources 



Source 


a 

(J2000.0) 


5 

(J2000.0) 


Flux Calibrator 


Gain Calibrator 


Distance 
(kpc) 


Mass 
(Af©) 


'ULSR 
(km s-i) 


G19.61-0.23'^ 


18 27 38.1 


-11 56 39.0 


Neptune 


1911-201 


3.50 


450 


40.0 


G29.96-0.02'' 


18 46 04.0 


-02 39 21.5 


Neptune 


1751-1-096 


6.0 


1100 


98.8 


IRAS 16293-2422<='' 


* 16 32 22.8 


-24 28 33.0 


MWC349 


1625-254 


0.16 


5.4 


3.9 










1733-130 









HFuruva et ahl l|2005h 
lOlmi et al] ll2003h 
" ISchoier et all l l2002h 
iDi Francesco et al.l ll200ll ') 



Table 2. OBSERVED MOLECULES AND TRANSITIONS 



Species 


Rotational Partition Function 


Transition 


Frequency (MHz) 


E^ (K) 


< 5j,jAj2 > (debye^) 


CHsCOOH^' 


Qr = U.lTr^^ 


10*,io -9*,9 E 


111,507.270(20) 


30.5 


54.8 




Qr = UAbTr^'^ 


10*,io -9,,9 A 


111,548.533(20) 


30.5 


54.8 


HCGOCHg'' 


92,8 ~ §2,7 E 


107,537.189(25) 


28.8 


22.8 






92,8 - 82,7 A 


107,543.746(25) 


28.8 


22.8 


GH3GH2CN<= 


12ll,2 — llll,l 


107,539.854 (3) 


167.8 


28.4 






12ii,i — llii,o 


107,539.854 (3) 


167.8 


28.4 






124,9 — 114,8 


107,543.918 (3) 


51.4 


158.1 






124,8 — 114,7 


107,547.593 (3) 


51.4 


158.1 






123,10 — 113,9 


107,594.040 (3) 


43.6 


166.8 



Note. — The \a uncertainty of the frequencies is in units of kHz. Each of the CH3COOH lines consists of two a- type and 
two b-type transitions. This is represented by an asterisk substituted for the K_ quantum numbers. 

Hllvushin et all ||2008|) 

lOesterling et all lll999l) 

' iFukuvama et al. I l ll996h : lRemiian et al.l ||2007| ') 
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Table 3. DETECTED MOLECULES AND TRANSITIONS TOWARD G19.61-0.23 



c 

op6C16S 


J_r3,iiSitioii 


AI (Jy beam ^) 


Av (km s "'") 


AT fill Knnm~l^ An 


^Km s J 






Fitted Data 


Model Data 




CH3COOH 


10^ in — 9* Q E 


0.07±0.01 


6.4±0.9 


0.10 


6.4 




10,, 10 - 9,,9 A 


0.09±0.02 


6.4±0.9 


0.12 


6.4 


HCOOCH3 


92,8 — 82,7 E 


0.20±0.02 


6.4±0.9 


0.18 


6.4 




92,8 - 82,7 A 


0.20±0.02 


6.4±0.9 


0.18 


6.4 


CH3CH2CN 


12ii,. - llii,. 


0.18±0.02 


12.9±0.4 


0.13 


12.9 




124,9 — 114,8 


0.59±0.02 


12.9±0.4 


0.66 


12.9 




124,8 — 114,7 


0.59±0.02 


12.9±0.4 


0.66 


12.9 




123,10 — 113,9 


0.71±0.01 


12.9±0.4 


0.72 


12.9 


U107.588 


107.5880 GHz 


0.16±0.04 


5.6±1.6 






U107.591 


107.5910 GHz 


0.27±0.03 


7.8±1.2 






U107.597 


107.5970 GHz 


0.31±0.03 


12.5±1.2 






U107.604 


107.6040 GHz 


0.23±0.03 


11.1±1.5 






U111.506 


111.5055 GHz 


0.06±0.02 


5.2±1.7 






U111.510 


111.5095 GHz 


0.07±0.01 


17.2it4.5 






U111.546 


111.5460 GHz 


0.07±0.01 


8.7±2.0 







Note. — The HCOOCH3 lines were fit assuming the same linewidth. The same line width fitting criteria 
applies for the CH3CH2CN and CH3COOH lines. For simplicity, the two CH3CH2CN fines listed in Table 
2 are now represented as one by an asterisk substituted for the K_ quantum numbers. The other three 
CH3CH2CN lines remain unchanged. 



Table 4. DETECTED MOLECULES AND TRANSITIONS TOWARD G29.96-0.02 



Species 


Transition 


AI (Jy beam-i) 


Av (km s~^) 


AI (Jy beam-i) 


Av (km s~^) 






Fitted Data 


Model Data 


CH3COOH 


10,, 10 - 9,, 9 E 


~0.04 










10,,io -9,,9 A 


<0.03 








HCOOCH3 


92,8 — 82,7 E 


0.08±0.01 


7.8±1.1 


0.09 


7.8 




92,8 - 82,7 A 


0.08±0.01 


7.8±1.1 


0.09 


7.8 


CH3CH2CN 


12ii,* — llii,* 


0.07±0.01 


10.8±0.4 


0.03 


10.8 




124,9 — 114,8 


0.21±0.01 


10.8±0.4 


0.23 


10.8 




124,8 — 114,7 


0.21±0.01 


10.8±0.4 


0.23 


10.8 




123,10 — 113,9 


0.27±0.01 


10.8±0.4 


0.26 


10.8 


U107.591 


107.5910 GHz 


O.llitO.Ol 


8.5±0.9 






U107.597 


107.5975 GHz 


0.08±0.01 


10.8±1.2 






U107.604 


107.6040 GHz 


O.lOitO.Ol 


12.2itl.3 






U111.541 


111.5410 GHz 


0.07±0.01 


12.7±2.0 







Note. — 



Table comments are the same as in Table 3. 
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Table 5. DETECTED MOLECULES AND TRANSITIONS TOWARD IRAS 16293-2422 



Species 


Transition 


AI (Jy beam-i) 


At) (km s -'-) 


AI (Jy beam~i) 


Ad 


(km s-i) 






Fitted Data 


Model Data 




IRAS 16293-2422 A 


CH3COOH 


10*, 10 - 9, ,9 E 


0.05±0.01 


6.0±2.0 


0.03 




6 




10,,io -9.,9 A 


0.05±0.01 


6.0±2.0 


0.05 




6 


HCOOCH3 


92,8-82,7 E 


0.050±0.009 


14.6±1.5 


0.049 




14.6 




92,8 — 82,7 A 


0.050±0.008 


14.6±1.5 


0.049 




14.6 


U111.508 


111.5080 GHz 


0.05±0.01 


5.0±2.0 








IRAS 16293-2422 B 


CH3COOH 


10*, 10 - 9, ,9 E 


<0.01 












10*,io -9*,9 A 


<0.01 










HCOOCH3 


92,8-82,7 E 


0.063±0.005 


5.5±0.5 


0.059 




5.5 




92,8 - 82,7 A 


0.063±0.005 


5.5±0.5 


0.059 




5.5 


U111.510 


111.5095 GHz 


0.06±0.01 


8.0±2.0 








U111.541 


111.5410 GHz 


0.06±0.01 


5.0±1.0 









Table 6. COLUMN DENSITIES OF THE OBSERVED MOLECULES 



Species (cm^^) G19.61-0.23 G29.96-0.02 I16293A I16293B 

CH3COOH (2.0±1.0)xl0i'^ < 9x10" ~1. 6x10^5 <6xl0" 

HCOOCH3 (9.0±2.0)xl0i6 (4.0±1.0)xl0i'' (1.5±0.3) x lO^^^ (6.0±1.0) x 10^^ 

CH3CH2CN (6.0±3.0)xl0i^ (I.l±0.7)xl0i6 <6xl0i=' <3xl0i* 



Table 7. COMPARISON OF CH3COOH SOURCES 



Species 


Sgr B2(N-LMH)=' 


W51e2'^ 


G34.3-f0.15= 


G19.61-0.23'^ 


IRAS 16293 A<i 


NCH3COOH (cm~^) 


-7.3x10^5 


1.7±0.5xl0^'^ 


(0.77-1.64) X 10^5 


2.0±1.0xl0^'' 


~1.6xl0^^ 


NCH3COOH/NHCOOCH3 


(4-7) X 10-2 


(1-6) X 10-2 


~3.3xl0-2 


2. 2±0. 1x10^1 


l.OitO.lxlO^i 



HMehringer et all lll997l) : Tr =200 K 
iRemiian eTsI\ hooit] : Tr =201 K 
teemiian et~ai] | |2003|) : Tr =70-185 K 

^this work; Tr =161 K for G19 and Tr =62 K for I16293A 
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(a) CH3COOH 




G19.61 — 0.23 


33" 








36" 








S (2000) 
K 






r 


42" 








ir56'45" 









(b) HCOOCH3 




o 



ia''27'"38;6 38?4 38f2 SSfO 37f8 

a (2000) 



18''27'"38?6 38?4 38f2 38;o 37;8 

a (2000) 



(c) CHjCHjCN 




ia''27"'38f6 



3a?4 



3B?2 
« (2000) 



3B?0 



37!8 



Fig. 1. — G19 maps, (a) CH3COOH contours overlaid on the gray continuum map. The continuum 
emission maps shown were made from channels contained in the wideband windows which were free 
from hne emission. The contour levels are -3, 3, 4 and 5 x d = 24 mJy beam^^. The gray scale 
unit is Jy beam~^. The synthesized beam size is l'.'65x 1'.'42 (P.A.=46.3°), shown in the bottom 
left corner. The continuum beam size is 1'.'58 x 1'.'45 (P.A.=36.4°). In the continuum map, there 
are two main components, A and C. The CH3COOH emission is mainly overlaid on component C. 
(b) HCOOCH3 contours overlaid on the gray continuum map. The contour levels averaged over 3 
channels arc 3, 6, 9, 12, 15, 18, 21, 24, 27 and 30 x cr = 20 mJy beam'^ The synthesized beam 
size is 2'.'16x 1'.'67 (P.A.=59.6°). (c) CH3CH2CN contours overlaid on the gray continuum map. 
The contour levels averaged over 6 channels, are 3, 9, 15, 21, 27, 33, 39, 45, 51, 57 and 63 x cr = 18 
mJy beam-^ The synthesized beam size is 2'.'16xl?67 (P.A.=60.3°). 
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G19.61 - 0.23 




(b) 














U111.506 U111.510 






Frequency (GHz) 



Fig. 2. — The observed spectra of G19.61-0.23. The first column (a) presents two HCOOCH3 lines, 
92,8 - 82,7 E and 92,8 - 82,7 A; and four CH3CH2CN lines, 124,9 - ll4,8, 124,8 - ll4,7, 12ii,2 - llii,i 
and 12ii^i — llii,o- The 12ii^2 — llii,i and 12ii^i — lln^o lines are not spectrally resolved. The v^sR 
range is -15 to 69.4 km s~^. The second column (b) contains only one known line, CH3CH2CN 
123,10 — 113,9- The vlsr range is -3.4 to 81 km s~^. The third column (c) contains the CH3COOH 
10*, 10 — 9*^9 E line and the fourth column (d) contains the CH3COOH 10*^io — 8*^9 A line. The rest 
frequencies of the windows are 107.54375, 107.59404, 111.50727 and 111.54853 GHz, respectively. 
The vlsr ranges are 11.3 to 92.7 and -1.8 to 79.6 km s^^ for column (c) and (d), respectively. The 
first row presents the observed spectra (Hanning smoothed over 3 channels) overlaid on modeled 
spectral line data represented by Gaussians. The second row shows the modeled spectral line 
data from the first row convolved with the spectral resolution of the original observation and with 
random noise added to the spectrum. The value of the random noise is obtained with the line-free 
continuum spectrum in the windows. Finally, the third row shows the residuals from the modeled 
data in the second row subtracted from the observational data in the first row. 
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a (2000) 

Fig. 3. — G29 maps, (a) HCOOCH3 contours overlaid on the gray continuum map. The contour 
levels averaged over 3 channels are -3, 3, 4, 5, 6, 7 and 8 x o" = 21 mJy beam^^. The synthesized 
beam size is 2'.'81 x 1'.'56 (P.A.=53.2°), shown in the bottom left corner, (b) CH3CH2CN contours 
overlaid on the gray continuum map. The contour levels averaged over 6 channels, are -3, 3, 5, 7, 
9, 11, 13 and 15 x cr = 16 mJy beam-^ The synthesized beam size is 2^'82xl?56 (P.A.=53.5°). 
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Fig. 4. — The observed spectra of G29.96-0.02. The arrangement of frequency and windows is the 
same as Figure El The ylsr ranges are 44.6 to 129, 56.2 to 140.6, 68.8 to 150.2 and 55.7 to 137.1 
km s~^, respectively. 
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16''32"'23?4 23?2 23?0 22?8 22?6 22?4 22?2 
a (2000) 

Fig. 5. — The HCOOCH3 contour maps are made individually for I16293A and B (denoted as A and 
B on the maps) due to the different linewidths. The contours are overlaid on the gray continuum 
map of IRAS 16293-2422. (a) The contours located at I16293A are made with the 92,8 - 82,7 E line 
averaged over six channels. The contour levels are -2, 2, 3, 4 and 5 x a = 21 mJy beam~^. (b) 
The contours located at I16293B are made with one channel of the 92,8 ~ 82,7 E line. The contour 
levels are -2, 2, 3, 4 and 5 x a = 25 mJy beam~^. 
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IRAS 16293-2422 A 




107.54 107.55 107.56 107.58 107. 5S 107.8 107.81 111.4S 111.5 111.51 111.53 111,54 111.55 111.58 



Frequency (GHz) 

Fig. 6. — The observed spectra of I16293A. The arrangement of frequency and windows is the same 
as FigureEl The column vlsr ranges are 44.6 to 129, 56.2 to 140.6, 68.8 to 150.2 and 55.7 to 137.1 
km s^^, respectively. 



IRAS 16293-2422 B 




107.54 107.55 107.56 107.58 107.59 107.6 107.61 111.49 111.5 111.51 111.53 111.54 111.55 111.56 

Frequency (GHz) 



Fig. 7. — The observed spectra of I16293B. The arrangement of frequency and windows is the same 
as Figure m The column vlsr ranges are -54.4 to 30, -42.7 to 41.6, -24.8 to 56.6 and -36.8 to 44.6 
km s~^, respectively. 
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Fig. 8.— IRAM 30-m observations toward G19.61-0.23 and G29.96-0.02 (|Fontani et al 
acetic acid lines are covered in this window at 111.50727 and 111.54853 GHz. 



20071). The 



